Dynamics of a three terminal mechanically flexible tunneling contact 
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The dynamics of a nanoelectromechanical system in the form of a three-terminal tunneling device 
is studied by analytical and numerical methods. The main results are the existence of bistable 
stationary states resulting in directly detectable chaotic behavior. 
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I. INTRODUCTION L 



The growing interest in nano electromechanical sys- 
tems (NEMS) is partly due to the advances in microfab- 
rication by self-assembly of biometallic composites where 
metallic or semiconducting clusters are combined with or- 
ganic molecules such as polymers or DNAErO. Since the 
organic chain molecules used in self-assembly are typi- 
cally a few orders of magnitude softermj than ordinary 
solids their utilization in nanoelectronics implies that me- 
chanical degrees of freedom can have great impact on the 
electronic transport properties in such systems. 

We have previouslyuli!3 investigated how the mechani- 
cal degrees of freedom couple to the electrical transport 
properties in a Coulomb blockade double junction where 
the central island is free to move in a parabolic poten- 
tial. A dynamical instability was shown to exist in this 
structure that caused the central island to oscillate be- 
tween the external electrodes. In the low temperature 
limit, where charging effects are important, this insta- 
bility gave rise to a mechanically mediated current that 
approached the value / = 2eNf where e is the elemen- 
tary charge, N the maximum number of excess electrons 
allowed on the island and / the frequency of elastic vibra- 
tions. The fluctuations of this current have been shown to. 
disappear exponentially with decreasing temperaturestil 
implying that this system may be well suited for current 
standard purposes. In realizing a current standard de- 
vice of this type a grain placed on a flexible cantilever 
positioned between the two electrodes can rexilacc the 
grain and the soft molecular links considered inBEj. Such 
nanoscale jnechanical resonators have successfully been 
fabricatedtj'Ej. It has also been suggested by Tuomi- 



nen et alS3, who studied a macroscopic electromechani- 
cal system at room temperature, that-carbon nanotubes 
connected to nanoscale metallic grainst3't9, could be used 
for this purpose. 

In this paper it is shown that when a system consisting 
of a grain attached to the tip of a cantilever by a tunnel 
junction (see Fig. |l]a) is treated as a three terminal device 
it shows a rich dynamical structure ranging from stable 
limit cycle behavior to deterministic chaos. Since this 
investigation focuses on the dynamical properties, the 
temperature is considered to be large enough for charging 
effects to be absent while still low enough for thermal 
fluctuations to be negligible. 
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FIG. 1. Schematic layout of the system, (a) A cantilever 
attached to a conducting grain through a tunnel junction sit- 
uated asymmetrically between two metallic electrodes sepa- 
rated by a distance L. The electrodes and the cantilever are 
biased with voltages Vl, Vr and Vb respectively. The dis- 
tance between the grain and the electrodes are such that elec- 
trons may tunnel through the structure, (b) Self assembled 
single electron transistor structure. The dynamical features 
discussed for system (a) are relevant also for this system. 

The same model also applies to the dynamics of the sin- 
gle electron transistor structure in Fig. [Tf3. The central 
island is here resting on a mechanically soft insulating 
monolayer allowing for center of mass motion of the grain. 
This system should of course be studied in the low tem- 
perature limit where charging effects are of importance. 
However, since much of the dynamical features of the sys- 
tem in the high temperature limit are expected to survive 
at lower temperatures the results obtained in this work 
is of interest for further studies in the low temperature 
regime. 
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II. MODEL SYSTEM 

The system, as depicted in Fig. |l|a consists of a small 
metallic grain, typically a few nanometers in diameter, 
attached to the tip of a metallic cantilever through a tun- 
nel junction. The other end of the cantilever is assumed 
to be clamped and connected to a voltage source Vb- 
The grain attached to the cantilever is situated between 
two metallic electrodes located close enough to the grain 
to allow tunneling but with tunneling resistances much 
larger than that of the grain-cantilever junction. If the 
typical frequency of the dynamical motion of the system 
is well below the plasmafrequency of the metallic compo- 
nents the charges and potentials on the conductors are 
related by the matrix Cy ; 
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FIG. 2. Electric circuit model. Each tunnel junction is 
modeled by a capacitor in parallel with a tunnel resistance 



Rl,r,b = G L 



Furthermore the grain-ground capacitance 



has been taken into account through the capacitance Cgo- 
Only the "nearest neighbor" capacitances are considered in 
this model. 



dt 



II + In 



If 



Using I l . r .b = (Vl.r.b - Vg)G l ,r,b where Glpjj are 
conductances of the tunnel junctions (cnf. Fig. |2j) one 
finds 



dQc 
dt 



= V L G L + V R G R + V B G B - V G (G L + Gr + G b ). 



The flexibility of the cantilever allows for a mechanical 
degree of freedom in the system. We will assume that this 
freedom is one-dimensional, i.e. that the cantilever only 
bends in one direction and that this bending is so small 
that we can consider the grain to move on a straight line 
between the two electrodes. The deflection 6 of the grain 
from the equilibrium position towards the right electrode 
is described by the equation of motion 

6 + 7 j + ujqS = £Qg I m, 

where loq is the elastic frequency of the system and 7 a 
parameter describing the damping. The mass m appear- 
ing on the right hand side is an effective mass depend- 
ing on the precise geometry and design of the cantilever- 
grain system. In the force term only the effect of the 
electro-static field between the left and right electrodes 
£ w (Vl — Vr)/L to linear order in Qq has been consid- 
ered. The tunneling conductances Gl.r have a sensitive 
dependence on the grain displacement 6 



Glm = G L R exp ( T- 



where the tunneling length A is determined by the work 
function cj) of the electrodes, 



2y / 2m e <fi 



13*3 ■ 



The system can be mapped onto an electrical network 
model seen in Fig. g. In this model only "nearest neigh- 
bor" coefficients in Cy have been taken into account. 
Each tunnel junction is modeled by a resistance R in 
parallel with a capacitance C the latter being related 
to the coefficeints Cy by linear transformations. Using 
this model, the grain potential Vq can be expressed in 
terms of the grain charge Qg and the applied bias volt- 
ages Vr, l ,b (cnf. Fig. 0) as 



V G (Qg) 



Qg + VrCr + V l C l + V b C b 



where Cs = Cl + Cr + Cb + Cgo- The currents Il,r,b 
flowing from the electrodes to the grain, as shown in 
Fig. ||, determine an equation of motion for Qg', 



The variations of the capacitances with position have 
been neglected since they are much smaller than the 
changes in the conductances which dominate the non- 
linear behavior of the system. By measuring the deflec- 
tion 6 in units of lambda i.e. S — £A, one arrives at the 
following system of equations, 



n = - 7 n- w &- 



Vl - Vr 
LmX 



Qc 



(i) 

(2) 



Q G = V L G° L e-i + VrG^ + V B G B - U G (Q G )G E (3) 
where 

G*{t) = Gla-t + G R et + G B (4) 

and 

T/ , n v Qg + VlC l + VrCr + V b Cb 

Vg(Qg) = ■ (5) 



2 




5 



FIG. 3. Solutions to the fixed point equation a£ = H(^). 
When ■& — Vl/Vb < weither one or three solutions exist 
depending on whether $ is larger or smaller than the ratio 
—Cl/Cb and what value a oc luq assumes. In (a) only one 
stationary point exists while in (b) either one or three sta- 
tionary points appear. 
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FIG. 4. Stationary operation. When the damping constant 
7 is large the static deflection of the cantilever will cause 
a pronounced transistor-like action due to the exponential 
decrease of the tunneling resistance between the grain and 
the right electrode. Figure (a) shows the current flowing from 
the grain to the right electrode as a function of the voltage Vl 
applied to the left electrode. The different curves correspond 
to different biases Vb applied to the bottom electrode. In (b) 
the displacement 5 of the cantilever and the current flowing 
from the left electrode to the grain is plotted for the case 
V B = 4.0 V. 



III. FIXED POINTS AND STABILITY 

In order to characterize the dynamical behavior of the 
system the existence of fixed points is investigated and 
then the stability of these points is considered. Through- 
out the rest of the article the right electrode potential will 
be used as reference i.e. we will put Vr = from here 
on. With this convention the stationary points satisfy 
the system of equations 

o = n 

o = _ 7 n- w gc + (^)Q G 

= V L G° L e-t + V B G B - Qg+Vl ^ +VbCb Gx(£)- 

(6) 

Defining a = uj 2 Lm\/C-£, r/ 2 — V% and er = VlVb, equa- 
tion (ra) can be recast in the form a£ = H(£) where 



r] 2 g L e * + a 
g L e~^ + gn^ + 1 



- (?? 2 CL + CTCb). 



Here the dimensionless conductances gu.L = G rl /Gb 
and capacitances cl,b — Cl,b/Cy. have been introduced. 
Since H(£) is a restricted function the equation = 
H(£) has at least one solution and at the most three 
solutions. Defining d = rf/a = Vl/Vb three different 
cases can be identified, 

• $ < 0, H(£) has one minimum located at £o = 
iln([l — $]<?l/3,r) which means that the system 
has either one or three fixed points. 

• < $ < 1, H(£) has one maximum located at 
£o = 7jln([l — , &\gL/gR) and again we may have 
either one or three fixed points. 

• $ > 1, In this case H(£) is monotonic hence only 
one fixed point can exist. 

When ■& < two different cases can be distinguished; 
i") < —cb/cl and > $ > —cb/cl- In the first case 
lim^^ +oc H(£) < and only one solution lying in the 
left half plane is possible. In the second case the corre- 
sponding limit is positive and if a is chosen small enough 
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three solution will appear, one in the left half plane and 
two in the right. The different scenarios are shown in 
Fig. ||. When < & < 1 one can again single out two 
cases; Cb/(1 — cl) < $ < 1 and < ?9 < cs/(l — cl). 
When cb/(1 - c L ) < ■& < 1 then lim^-oo H(£) > and 
only one solution located in the right half plane is possi- 
ble. In the other case the limit is negative and one can 
find two more solutions in the left half plane by choosing 
a sufficiently small. 

A stability analysis of the fixed points obtained above 
shows that: In the case of only one fixed point this point 
will be a stable node, i.e. all-eigenvalues Jacobian of the 
system are real and negativaLU, if the damping 7 exceeds 
the critical damping 7 C , 

y(r-i(0+c 2 r(0) 2 -4^^) • (7) 
Here the total i?C-timc t(£) is defined as 

In the case of three fixed points one will be condition- 
ally stable depending on whether 7 is larger or smaller 
than 7 C . The second "middle" one (cnf. Fig. |b) will 
not be a node but instead be a saddle point of index 1 
(two negative and one positive eigenvalue of the Jaco- 
bian) and hence not a stationary point irrespectively of 
the value of 7 while the remaining point (corresponding 
to the rightmost solution in figure gp) will always be a 
stable node. 



R° = 1e+09O_ 
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FIG. 5. Bistable operation. When a is small enough two 
stable fixed points emerge leading to a bistable situation. The 
solid line corresponds to the "expected" solution while the 
dashed lines correspond to the new second stable root to the 
fixed-point equations. 



IV. STATIONARY OPERATION 

When y^a, which is proportional to the square of the 
frequency of elastic vibrations ojq, is large compared to 
the applied bias voltages one expects to find only one 
stationary solution as discussed above. In Fig. [|a the I — 
V characteristics for this case is shown for an asymmetric 
setup with R° R = 1 Gft, R° L = 1 TCI, R B = 6 MQ and 
all capacitances set to 1 aF. The figure shows the current 
flowing from the grain to the right electrode (as opposed 
to Ir defined in Fig. |^ which was defined in the opposite 
direction) as a function of the voltage applied to the left 
electrode. The different curves correspond to different 
biases applied to the bottom electrode. For negative Vl 
almost no current flows in the system since the grain is 
essentially disconnected from both leads. As the field is 
increased the cantilever will start to deflect towards the 
right electrode causing an exponential decrease in Rl 
allowing the current to grow. 

As Vl is further increased the charge on the grain will 
eventually become negative due to the capacitive cou- 
pling to the left electrode resulting in a decreased deflec- 
tion disconnecting it from the leads once again. In Fig. |] 
the the displacement of the grain along with the current 
flowing to it from the left electrode is shown as a func- 
tion of Vl when Vb — 4.0 V. This current is essentially 
zero until a bias of approximately —3.0 V is reached. 
Furthermore, from this graph it can be seen that the dis- 
placement of the grain is just a few times A (for Au A is 
typically 0.5 A). 

For i/a small compared to the applied bias voltages, it 
is possible to have two stable fixed points in the system. 
Using exactly the same parameters as above but with a 
reduced a this second solution appears. In Fig. || the 
I —V characteristics of this bistable mode is shown. The 
solid line correspond to the same solution as above and 
the dashed solution to the new one that appear due to 
the reduced a. 
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FIG. 6. Phase space trajectory for a stable limit cycle. As 
the damping 7 is reduced below the critical damping 7 C the 
system settles in to a stable limit cycle. 
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V. DYNAMICAL OPERATON 



A. One Fixed Point 



In the dynamical regime the system displays a rich 
structure. One of the most interesting features is that 
this is a nanoscale system with directly detectable chaotic 
behavior. This means that in order to determine what 
type of motion the system exhibits it is sufficient to mon- 
itor the currents. Due to the multitude of parameters 
and the system's complex dependence on these, only a 
few archetypical cases will be illustrated by means of nu- 
merical integration of the equations of motion (|l|)-(^). In 
the presented simulations a system with the same pa- 
rameters as in the static case (see the previous section) 
is considered but with a damping rate 7 reduced below 
7 C . The different simulations then correspond to different 
sets of bias voltages Vb and Vl ■ 
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FIG. 7. Multiply period doubled phase-space trajectory 
and the Corresponding plot of the currents II and Ib- (a) Bi- 
asing the system towards the situation with three fixed points 
will cause subsequent period doublings of the limit cycle, (b) 
By monitoring the currents flowing in the left lead and the 
bottom lead at the same time these period doublings can be 
detected. 



We first consider the situation when the system has 
only one fixed point corresponding to the situation in 
Fig. ||a. This is achieved by using a fixed voltage Vb = 3 
V and imposing a positive bias voltage Vl ■ For small val- 
ues of this voltage the system remains stable, as expected, 
until a critical threshold voltage is reached. Further bias- 
ing leads, for positive Vl, to a limit cycle regime. For the 
range of positive voltages where the algorithm was sta- 
ble this cycle remained. For negative bias voltages there 
exists a threshold voltage as well, and as this is reached, 
a stable limit cycle appears. An example of this type of 
motion is shown in Fig. || (recorded at Vl — —0.35 V). 
Decreasing Vl moves the system toward the situation in 
Fig. |^b i.e. we approach the situation with three fixed 
points. This leads to a sequence of period doublings. The 
behavior in this regime is illustrated in Fig. ^a. These pe- 
riod doublings can be directly detected by simultaneously 
measuring the currents Is and II and plotting them as 
in Fig. 0b. Lowering the bias more eventually leads to a 
totally chaotic regime like the one in Fig. || (Vl = —0.6 
V) which is again reflected in the currents. 
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FIG. 8. Chaotic motion, (a) Phase space trajectory in the 
chaotic regime. Biasing the system very close to the situation 
with three fixed points the system becomes chaotic, (b) Chaos 
is also reflected in the corresponding plot of the currents II 
and Ib- 



5 



Further lowering of Vl after this point leads to an alter- 
nating series of period doubled limit cycles and chaotic 
trajectories until three fixed points appear in the system. 



B. Three Fixed Points 

In order to be able to study the case with three fixed 
points, Vb was raised to 6 V while Vl was set to 0.7 
V. This corresponds to the situation with three solutions 
in Fig. |^b. Numerical integration revealed the structure 
displayed in Fig. ||. Starting close to the conditionally 
stable fixed point, the leftmost one in Fig. ||b, a stable 
limit cycle is eventually reached. Starting the simulation 
in the vicinity of the middle one (always unstable) the 
trajectory cither connects to this limit cycle or becomes 
attracted by the third fixed point, which is always stable. 
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FIG. 9. Phase space trajectory when three fixed points are 
present. When the system has three fixed points, one of them 
will always be stable, one will always be unstable and one will 
be is unstable if 7 < j c . The three fixed points are indicated 
by * in the figure. The stable limit cycle that exists in this 
case can also be seen in the figure. 



VI. CONCLUSIONS 

We have shown that the three terminal flexible tun- 
neling structures in Fig. [j], which are of interest for both 
current standard purposes as for self assembled quantum 
devices, have several characteristic dynamical features: 
When the damping 7 in the system is high (low quality 
factor) the system displays a stationary behavior which 
for some parameter values can be bistable. If the quality 
factor is large enough, i.e. the damping satisfies 7 < j c , 
the dynamics of the system range from stable limit cycle 
behavior to deterministic chaos. It is furthermore possi- 
ble to have a situation where stable fixed points coexist 
with stable limit cycles. The chaotic motion of the sys- 
tem can be directly detected by measuring the currents 
flowing from the terminals. 
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